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We propose a new non-perturbative method to compute derivatives of gauge coupling constants with respect to 
anisotropic lattice spacings (anisotropy coefficients). Our method is based on a precise measurement of the finite 
temperature deconfining transition curve in the lattice coupling parameter space extended to anisotropic lattices 
by applying the spectral density method. We determine the anisotropy coefficients for the cases of SU{2) and 
SU (3) gauge theories. A longstanding problem, when one uses the perturbative anisotropy coefficients, is a non¬ 
vanishing pressure gap at the deconfining transition point in the SU{3) gauge theory. Using our non-perturbative 
anisotropy coefficients, we hud that this problem is completely resolved. 


1. Introduction 

In a phenomenological study of heavy ion col¬ 
lisions and evolution of early Universe, it is im¬ 
portant to evaluate the energy density and the 
pressure of the quark-gluon plasma near the tran¬ 
sition temperature of the deconfining phase tran¬ 
sition of QCD. 

On an anisotropic lattice with Os and at the lat¬ 
tice spacings in spatial and temporal directions, 
the standard plaquette action for SU{Nc) gauge 
theory is given by S' = -/J* X] Ps{x)-ht L Pt{.x), 
where Ps{t) is the spatial (temporal) plaquette. 
Hence the energy density, e = — -h , and the 
pressure, p = are expressed in terms of 

the anisotropy coefficients, 

We choose at and ^ = as/at as independent vari¬ 
ables to vary the lattice spacings. 

Perturbative values for these anisotropy coef¬ 
ficients were calculated by Karsch in Q|. How¬ 
ever, when we apply them to data obtained 
by MC simulations, we encounter pathologi¬ 
cal results such as a negative pressure and a 
non-vanishing pressure gap at the deconfining 
transition point of SU{3) gauge theory. Non- 
perturbative anisotropy coefficients are, there¬ 
fore, required to study e and p in MC simulations. 

Two non-perturbative methods have been 
adopted to determine the anisotropy coefficients. 
One is “the matching method” based on a 


measurement of ^ as a function of Ps and /?* by 
matching spatial and temporal Wilson loops. The 
other is a method based on a non-perturbative 
estimate of pressure obtained by “the integral 
method” [|§. 

In this paper, we propose a new non- 
perturbative approach to compute the anisotropy 
coefficients, and determine the coefficients for 
SU{2) and SU{3) gauge theories [Q. We re¬ 
strict ourselves to the case of isotropic lattices, 
Ps = Pt = P, where most simulations are per¬ 
formed. In this case, two anisotropy coefficients 
are just the beta-function at 5 = 1; (at^|j)j=i = 
{at^^)^=i = whose non-perturbative val¬ 

ues are well studied both in SU{2) and SU{3) 
gauge theories . Furthermore, a combina¬ 

tion of the remaining two anisotropy coefficients is 
known to be again related to the beta-function by 

i- Therefore, only 

one additional input is required to determine the 
anisotropy coefficients for isotropic lattices. 


2. Method 

Our method is based on an observation that, 
the transition temperature Tc = 1/{Ntat{Ps, Pt)} 
is independent of the anisotropy of the lattice. 
This brings us the following relation between the 
anisotropy coefficients and the slope rt of the 
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transition curve in the {f3s, Pt) plane at ^ = 1, 


d/3, {dd^=J\dd^=i 

From this equation, we obtain the expressions 
for the customarily used forms for the anisotropy 
coefficients Cs(t) = {dg~^t^/dOas-.fi^ed where 
Ps = and /3, = 


Cs 


Ct 


1 fa n-2 

2Nc 2(1 + r,) °da J ’ 

1 f . , 1 - 2r, ^ 

2Nc\ 2{l + rt)^da 


( 2 ) 


Therefore, when the value for the beta-function 
is available, we can determine these anisotropy 
coefficients by measuring r, from the finite tem¬ 
perature transition curve in the {Ps, Pt) plane. 

As a result, e and p are given by 


e-3p 4 d/3 

—— = — SA/. a-—j 
r4 * da ‘ 


-p 




+ (Pt) - 2(P)o}, (3) 

= 3iV4a^^{(P,)-(P,)}, (4) 

da r, -I- 1 


where {P)o is the plaquette at T = 0. 

In order to determine the transition curve 
in the {Ps,Pt) plane, we compute the rotated 
Polyakov loop L. We define the transition point 
as the peak position of the susceptibility x = 
N^({L'^) — (T)^). The coupling parameter depen¬ 
dence of X in the (Ps,Pt) plane is computed by ap¬ 
plying the spectral density method |j^ extended 
to anisotropic lattices. This enables us to com¬ 
pute the anisotropy coefficients directly from sim¬ 
ulations at ^ w 1 without introducing an interpo¬ 
lation Ansatz, unlike the case of previous studies. 


3. Results 

We first test the method for the case of SU (2) 
gauge theory at the transition point Pc for W = 4 
and 5. Simulations are performed on 16^ x 4 and 
20 ^ X 5 lattices. Results for Cg and c, are de¬ 
noted by filled circles in Fig. (top). Our results 
are consistent with the results from the integral 
method (doted curves) ||]. 

We then study the more realistic case of the 
SU{3) gauge theory. Because the method works 




Figure 1. Anisotropy coefficients Cg and c, for 
the SU(2) and SU(3) gauge theories. 


well even with data obtained only on isotropic 
lattices, we analyze the high statistics data by the 
QCDPAX Collaboration ||^. Simulations were 
performed at the deconfining transition point for 
Nt = 4 and 6 on five lattices. Details of the SU (3) 
simulations are given in [p4| . 

Fig. § shows the (/3s, Pt) dependence of the sus¬ 
ceptibility on a 24^ x 36 x 4 lattice. Because 
the peak of the susceptibility becomes sharper as 
the spatial volume of the lattice is increased, we 
can measure r, most precisely on the spatially 
largest lattices. Therefore, in the following, we 
use the results obtained on the largest 24^ x 36 x 4 
and 36^ x 48 x 6 lattices. The values obtained 
on smaller lattices are consistent. For the beta- 
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Figure 2. The Polyakov loop susceptibility for 
SU{3) obtained on a 24^ x 36 x 4 lattice. 


Table 1 

Latent heat, pressure gap, A(Pt) /A(Ps) and — 
at the deconfining transition point of SU (3) gauge 
theory. 


lattice 

f3 

24^ X 36 X 4 
5.6925 

36^ X 48 X 6 
5.8936 

Ae/T^ 

2.074(34) 

1.569(40) 

Ap/p4 

0.001(15) 

-0.003(17) 

A(P} /A(P.) 

1.201(35) 

1.218(46) 

-n 

1.201(1) 

1.220(3) 


function, we adopt a result computed from a re¬ 
cent string tension data |^. 

In Fig. 1^ (bottom), we summarize our results 
for the Cs and ct of the SU (3) gauge theory (filled 
circles) together with previous values: the pertur¬ 
bative results (dot-dashed lines) |^, results from 
the integral method (doted curves) , and those 
from the matching of Wilson loops on anisotropic 
lattices (squares Q, triangles 0). We find that 
all non-perturbative methods give values which 
are roughly consistent with each other, showing a 
clear deviation from the perturbation theory. 

The deconfining transition is of first order for 
SU{S). At a first order transition point, we have 
a finite gap for energy density, the latent heat, 
but expect no gap for pressure. It is known 
that the perturbative anisotropy coefficients leads 
to a non-vanishing pressure gap at the decon¬ 


fining transition point: Ap/T'^ = —0.32(3) and 
-0.14(2) at At = 4 and 6 0. 

New values for the gaps in e and p using our 
non-perturbative anisotropy coefficients are sum¬ 
marized in Table ^ We find that the problem of 
non-zero pressure gap is completely resolved with 
our non-perturbative anisotropy coefficients. 

We note that, because the beta-function ap¬ 
pears only as a common overall factor in (^) and 
the conclusion that Ap vanishes with our 
anisotropy coefficients does not depend on the 
value of the beta-function. Actually we have from 
eqs.(^ and (^) a simple condition for Ap = 0: 

A{Pt) /A(P,) = -rt (5) 

where A(Ps(t)) is the gap in the spatial (tempo¬ 
ral) plaquette between the two phases. Although 
the two sides of (H) are obtained from quite differ¬ 
ent measurements, they agree precisely with each 
other as shown in Table |^. 
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